Large-Scale RT-qPCR Diagnostics for Seed Potato Certification.

Every year, Agroscope examines nearly 300,000 tubers for the presence of viruses, as regulated for the certification of seed potatoes intended for Swiss growers. Since 2016, this examination has been performed via RT-qPCR on dormant tubers directly after harvest. This method offers fast results and eliminates the need for the use of Rindite, which is a toxic and polluting gaseous compound previously used in Switzerland to break the dormancy of seed tubers. The implementation of this molecular analytical method for the routine diagnosis of regulated viruses makes it possible to conduct additional analyses via Illumina sequencing to assess the conformity of the primers and probes used with the sequences of the different viral isolates. This form of quality control in routine diagnosis is a source of information that can answer more fundamental scientific questions related to the epidemiology of viral strains related to certification. The datasets produced in this framework can also be used to explore the diversity of rare or unknown virus species in potato crops. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11540-021-09491-3.

Biotic and abiotic factors shape arbuscular mycorrhizal fungal communities associated with the roots of the widespread fern Botrychium lunaria (Ophioglossaceae).

Arbuscular mycorrhizal fungi (AMF) play central roles in terrestrial ecosystems by interacting with both above and belowground communities as well as by influencing edaphic properties. The AMF communities associated with the roots of the fern Botrychium lunaria (Ophioglossaceae) were sampled in four transects at 2400 m a.s.l. in the Swiss Alps and analyzed using metabarcoding. Members of five Glomeromycota genera were identified across the 71 samples. Our analyses revealed the existence of a core microbiome composed of four abundant Glomus operational taxonomic units (OTUs), as well as a low OTU turnover between samples. The AMF communities were not spatially structured, which contrasts with most studies on AMF associated with angiosperms. pH, microbial connectivity and humus cover significantly shaped AMF beta diversity but only explained a minor fraction of variation in beta diversity. AMF OTUs associations were found to be significant by both cohesion and co-occurrence analyses, suggesting a role for fungus-fungus interactions in AMF community assembly. In particular, OTU co-occurrences were more frequent between different genera than among the same genus, rising the hypothesis of functional complementarity among the AMF associated to B. lunaria. Altogether, our results provide new insights into the ecology of fern symbionts in alpine grasslands.

Small RNA-Omics for Virome Reconstruction and Antiviral Defense Characterization in Mixed Infections of Cultivated Solanum Plants.

In plants, RNA silencing-based antiviral defense generates viral small RNAs (sRNAs) faithfully representing the viral genomes. We employed sRNA sequencing and bioinformatics (sRNA-omics) to characterize antiviral defense and to reconstruct the full genomic sequences and their variants in the evolving viral quasispecies in cultivated solanaceous plants carrying mixed infections. In naturally infected Solanum tuberosum (potato), one case study revealed a virome comprising Potato virus Y (genus Potyvirus) and Potato virus X (genus Potexvirus), which was reconstructed by de novo-assembling separate genome-size sRNA contigs. Another case study revealed a virome comprising NTN and O strains of Potato virus Y, whose sRNAs assembled in chimeric contigs, which could be disentangled on the basis of reference genome sequences. Both viromes were stable in vegetative potato progeny. In a cross-protection trial of Solanum lycopersicum (tomato), the supposedly protective mild strain CH2 of Pepino mosaic virus (genus Potexvirus) was tested for protection against strain LP of the same virus. Reciprocal mechanical inoculations eventually resulted in co-infection of all individual plants with CH2 and LP strains, reconstructed as separate sRNA contigs. LP invasions into CH2-preinfected plants and vice versa were accompanied by alterations of consensus genome sequences in viral quasispecies, indicating a potential risk of cross-protection measures. Additionally, the study also revealed, by reconstruction from sRNAs, the presence of the mechanically nontransmissible Southern tomato virus (genus Amalgavirus) in some plants. Our in-depth analysis of sRNA sizes, 5'-nucleotide frequencies and hotspot maps revealed similarities in sRNA-generating mechanisms in potato and tomato, differential silencing responses to virome components and potential for sRNA-directed cross-targeting between viral strains which could not, however, prevent the formation of stable viromes.

Ultra-low input transcriptomics reveal the spore functional content and phylogenetic affiliations of poorly studied arbuscular mycorrhizal fungi.

Arbuscular mycorrhizal fungi (AMF) are a group of soil microorganisms that establish symbioses with the vast majority of land plants. To date, generation of AMF coding information has been limited to model genera that grow well axenically; Rhizoglomus and Gigaspora. Meanwhile, data on the functional gene repertoire of most AMF families is non-existent. Here, we provide primary large-scale transcriptome data from eight poorly studied AMF species (Acaulospora morrowiae, Diversispora versiforme, Scutellospora calospora, Racocetra castanea, Paraglomus brasilianum, Ambispora leptoticha, Claroideoglomus claroideum and Funneliformis mosseae) using ultra-low input ribonucleic acid (RNA)-seq approaches. Our analyses reveals that quiescent spores of many AMF species harbour a diverse functional diversity and solidify known evolutionary relationships within the group. Our findings demonstrate that RNA-seq data obtained from low-input RNA are reliable in comparison to conventional RNA-seq experiments. Thus, our methodology can potentially be used to deepen our understanding of fungal microbial function and phylogeny using minute amounts of RNA material.

Revisiting the Roles of Tobamovirus Replicase Complex Proteins in Viral Replication and Silencing Suppression.

Tobamoviral replicase possesses an RNA-dependent RNA polymerase (RDR) domain and is translated from genomic (g)RNA via a stop codon readthrough mechanism at a one-to-ten ratio relative to a shorter protein lacking the RDR domain. The two proteins share methyltransferase and helicase domains and form a heterodimer implicated in gRNA replication. The shorter protein is also implicated in suppressing RNA silencing-based antiviral defenses. Using a stop codon mutant of Oilseed rape mosaic tobamovirus (ORMV), we demonstrate that the readthrough replicase (p182) is sufficient for gRNA replication and for subgenomic RNA transcription during systemic infection in Nicotiana benthamiana and Arabidopsis thaliana. However, the mutant virus displays milder symptoms and does not interfere with HEN1-mediated methylation of viral short interfering (si)RNAs or plant small (s)RNAs. The mutant virus tends to revert the stop codon, thereby restoring expression of the shorter protein (p125), even in the absence of plant Dicer-like activities that generate viral siRNAs. Plant RDR activities that generate endogenous siRNA precursors do not prevent replication or movement of the mutant virus, and double-stranded precursors of viral siRNAs representing the entire virus genome are likely synthesized by p182. Transgenic expression of p125 partially recapitulates the ORMV disease symptoms associated with overaccumulation of plant sRNAs. Taken together, the readthrough replicase p182 is sufficient for viral replication and transcription but not for silencing suppression. By contrast, the shorter p125 protein suppresses silencing, provokes severe disease symptoms, causes overaccumulation of unmethylated viral and plant sRNAs but it is not an essential component of the viral replicase complex.

Low number of fixed somatic mutations in a long-lived oak tree.

Because plants do not possess a defined germline, deleterious somatic mutations can be passed to gametes, and a large number of cell divisions separating zygote from gamete formation may lead to many mutations in long-lived plants. We sequenced the genome of two terminal branches of a 234-year-old oak tree and found several fixed somatic single-nucleotide variants whose sequential appearance in the tree could be traced along nested sectors of younger branches. Our data suggest that stem cells of shoot meristems in trees are robustly protected from the accumulation of mutations.

Technical Insights into Highly Sensitive Isolation and Molecular Characterization of Fixed and Live Circulating Tumor Cells for Early Detection of Tumor Invasion.

Circulating Tumor Cells (CTC) and Circulating Tumor Microemboli (CTM) are Circulating Rare Cells (CRC) which herald tumor invasion and are expected to provide an opportunity to improve the management of cancer patients. An unsolved technical issue in the CTC field is how to obtain highly sensitive and unbiased collection of these fragile and heterogeneous cells, in both live and fixed form, for their molecular study when they are extremely rare, particularly at the beginning of the invasion process. We report on a new protocol to enrich from blood live CTC using ISET® (Isolation by SizE of Tumor/Trophoblastic Cells), an open system originally developed for marker-independent isolation of fixed tumor cells. We have assessed the impact of our new enrichment method on live tumor cells antigen expression, cytoskeleton structure, cell viability and ability to expand in culture. We have also explored the ISET® in vitro performance to collect intact fixed and live cancer cells by using spiking analyses with extremely low number of fluorescent cultured cells. We describe results consistently showing the feasibility of isolating fixed and live tumor cells with a Lower Limit of Detection (LLOD) of one cancer cell per 10 mL of blood and a sensitivity at LLOD ranging from 83 to 100%. This very high sensitivity threshold can be maintained when plasma is collected before tumor cells isolation. Finally, we have performed a comparative next generation sequencing (NGS) analysis of tumor cells before and after isolation from blood and culture. We established the feasibility of NGS analysis of single live and fixed tumor cells enriched from blood by our system. This study provides new protocols for detection and characterization of CTC collected from blood at the very early steps of tumor invasion.

Cell-free DNA testing of an extended range of chromosomal anomalies: clinical experience with 6,388 consecutive cases.

PURPOSE: Cell-free DNA (cfDNA) testing for fetal aneuploidies was broadly implemented for common trisomies and sex-chromosome anomalies (SCAs). However, such an approach identifies only 75 to 85% of clinically relevant aneuploidies. METHODS: We present a consecutive series of 6,388 cases, thus uncovering a broader array of aneuploidies, including the rare autosomal trisomies (RATs) and the maternally inherited deletion and duplication copy-number variations (CNVs), with complete and stratified follow-up by amniocentesis. Combined measurements of z-scores and the fetal fraction, in conjunction with fetal cfDNA enrichment, were used to stratify the likelihood of true and false results. RESULTS: We obtained an incremental diagnostic yield of 50%; RATs and CNVs were found to be significant causes of fetal pathology. Scrutinizing z-scores and the fetal fraction made it possible to distinguish the sources of false-negative results; predict the likelihood of false-positive results for major trisomies and SCAs; classify maternal mosaic SCAs and CNVs, preventing false-positive results; and robustly identify maternally inherited CNVs and detect recurrent genomic disorders as a standardized function of the fetal fraction. CONCLUSION: With the clinical pertinence of this broader detection scheme confirmed, we offer recommendations for its implementation.Genet Med 19 2, 169-175.

Roles of Non-Coding RNA in Sugarcane-Microbe Interaction.

Studies have highlighted the importance of non-coding RNA regulation in plant-microbe interaction. However, the roles of sugarcane microRNAs (miRNAs) in the regulation of disease responses have not been investigated. Firstly, we screened the sRNA transcriptome of sugarcane infected with Acidovorax avenae. Conserved and novel miRNAs were identified. Additionally, small interfering RNAs (siRNAs) were aligned to differentially expressed sequences from the sugarcane transcriptome. Interestingly, many siRNAs aligned to a transcript encoding a copper-transporter gene whose expression was induced in the presence of A. avenae, while the siRNAs were repressed in the presence of A. avenae. Moreover, a long intergenic non-coding RNA was identified as a potential target or decoy of miR408. To extend the bioinformatics analysis, we carried out independent inoculations and the expression patterns of six miRNAs were validated by quantitative reverse transcription-PCR (qRT-PCR). Among these miRNAs, miR408-a copper-microRNA-was downregulated. The cleavage of a putative miR408 target, a laccase, was confirmed by a modified 5'RACE (rapid amplification of cDNA ends) assay. MiR408 was also downregulated in samples infected with other pathogens, but it was upregulated in the presence of a beneficial diazotrophic bacteria. Our results suggest that regulation by miR408 is important in sugarcane sensing whether microorganisms are either pathogenic or beneficial, triggering specific miRNA-mediated regulatory mechanisms accordingly.

Microbial Communities of Conducting and Respiratory Zones of Lung-Transplanted Patients.

Background: Lung transplantation (LT) is a recognized treatment for end-stage pulmonary disease. Bacteria from the recipient nasopharynx seed the new lungs leading to infections and allograft damage. Understanding the characteristics and topological variations of the microbiota may be important to apprehend the pathophysiology of allograft dysfunction. Objectives: To examine the characteristics and relationship of bacterial compositions between conducting and respiratory zones of the allograft. Methods: We performed 16S rRNA gene sequencing on bronchial aspirates (BAs) and bronchoalveolar lavages (BALs) collected in pairs in 19 patients at several time-points post-LT. Results: The respiratory zone was characterized independently of the time post-LT by a higher bacterial richness than the conducting zone (p = 0.041). The phyla Firmicutes and Proteobacteria dominated both sampling zones, with an inverse correlation between these two phyla (Spearman r = -0.830). Samples of the same pair, as well as pairs from the same individual clustered together (Pseudo-F = 3.8652, p < 0.01). Microbiota of BA and BAL were more closely related in samples from the same patient than each sample type across different patients, with variation in community structure being mainly inter-individual (p < 0.01). Both number of antibiotics administered (p < 0.01) and time interval post-LT (p < 0.01) contributed to the variation in global microbiota structure. Longitudinal analysis of BA-BAL pairs of two patients showed dynamic wave like fluctuations of the microbiota. Conclusions: Our results show that post-transplant respiratory zones harbor higher bacterial richness, but overall similar bacterial profiles as compared to conductive zones. They further support an individual microbial signature following LT.

Evidence for the sexual origin of heterokaryosis in arbuscular mycorrhizal fungi.

Sexual reproduction is ubiquitous among eukaryotes, and fully asexual lineages are extremely rare. Prominent among ancient asexual lineages are the arbuscular mycorrhizal fungi (AMF), a group of plant symbionts with a multinucleate cytoplasm. Genomic divergence among co-existing nuclei was proposed to drive the evolutionary success of AMF in the absence of sex(1), but this hypothesis has been contradicted by recent genome analyses that failed to find significant genetic diversity within an AMF isolate(2,3). Here, we set out to resolve issues surrounding the genome organization and sexual potential of AMF by exploring the genomes of five isolates of Rhizophagus irregularis, a model AMF. We find that genetic diversity in this species varies among isolates and is structured in a homo-dikaryon-like manner usually linked with the existence of a sexual life cycle. We also identify a putative AMF mating-type locus, containing two genes with structural and evolutionary similarities with the mating-type locus of some Dikarya. Our analyses suggest that this locus may be multi-allelic and that AMF could be heterothallic and bipolar. These findings reconcile opposing views on the genome organization of these ubiquitous plant symbionts and open avenues for strain improvement and environmental application of these organisms.

The Bee Microbiome: Impact on Bee Health and Model for Evolution and Ecology of Host-Microbe Interactions.

As pollinators, bees are cornerstones for terrestrial ecosystem stability and key components in agricultural productivity. All animals, including bees, are associated with a diverse community of microbes, commonly referred to as the microbiome. The bee microbiome is likely to be a crucial factor affecting host health. However, with the exception of a few pathogens, the impacts of most members of the bee microbiome on host health are poorly understood. Further, the evolutionary and ecological forces that shape and change the microbiome are unclear. Here, we discuss recent progress in our understanding of the bee microbiome, and we present challenges associated with its investigation. We conclude that global coordination of research efforts is needed to fully understand the complex and highly dynamic nature of the interplay between the bee microbiome, its host, and the environment. High-throughput sequencing technologies are ideal for exploring complex biological systems, including host-microbe interactions. To maximize their value and to improve assessment of the factors affecting bee health, sequence data should be archived, curated, and analyzed in ways that promote the synthesis of different studies. To this end, the BeeBiome consortium aims to develop an online database which would provide reference sequences, archive metadata, and host analytical resources. The goal would be to support applied and fundamental research on bees and their associated microbes and to provide a collaborative framework for sharing primary data from different research programs, thus furthering our understanding of the bee microbiome and its impact on pollinator health.

Unusual RNA plant virus integration in the soybean genome leads to the production of small RNAs.

Horizontal gene transfer (HGT) is known to be a major force in genome evolution. The acquisition of genes from viruses by eukaryotic genomes is a well-studied example of HGT, including rare cases of non-retroviral RNA virus integration. The present study describes the integration of cucumber mosaic virus RNA-1 into soybean genome. After an initial metatranscriptomic analysis of small RNAs derived from soybean, the de novo assembly resulted a 3029-nt contig homologous to RNA-1. The integration of this sequence in the soybean genome was confirmed by DNA deep sequencing. The locus where the integration occurred harbors the full RNA-1 sequence followed by the partial sequence of an endogenous mRNA and another sequence of RNA-1 as an inverted repeat and allowing the formation of a hairpin structure. This region recombined into a retrotransposon located inside an exon of a soybean gene. The nucleotide similarity of the integrated sequence compared to other Cucumber mosaic virus sequences indicates that the integration event occurred recently. We described a rare event of non-retroviral RNA virus integration in soybean that leads to the production of a double-stranded RNA in a similar fashion to virus resistance RNAi plants.

MISIS-2: A bioinformatics tool for in-depth analysis of small RNAs and representation of consensus master genome in viral quasispecies.

In most eukaryotes, small RNA (sRNA) molecules such as miRNAs, siRNAs and piRNAs regulate gene expression and repress transposons and viruses. AGO/PIWI family proteins sort functional sRNAs based on size, 5'-nucleotide and other sequence features. In plants and some animals, viral sRNAs are extremely diverse and cover the entire viral genome sequences, which allows for de novo reconstruction of a complete viral genome by deep sequencing and bioinformatics analysis of viral sRNAs. Previously, we have developed a tool MISIS to view and analyze sRNA maps of viruses and cellular genome regions which spawn multiple sRNAs. Here we describe a new release of MISIS, MISIS-2, which enables to determine and visualize a consensus sequence and count sRNAs of any chosen sizes and 5'-terminal nucleotide identities. Furthermore we demonstrate the utility of MISIS-2 for identification of single nucleotide polymorphisms (SNPs) at each position of a reference sequence and reconstruction of a consensus master genome in evolving viral quasispecies. MISIS-2 is a Java standalone program. It is freely available along with the source code at the website http://www.fasteris.com/apps.

Field Trial and Molecular Characterization of RNAi-Transgenic Tomato Plants That Exhibit Resistance to Tomato Yellow Leaf Curl Geminivirus.

RNA interference (RNAi) is a widely used approach to generate virus-resistant transgenic crops. However, issues of agricultural importance like the long-term durability of RNAi-mediated resistance under field conditions and the potential side effects provoked in the plant by the stable RNAi expression remain poorly investigated. Here, we performed field trials and molecular characterization studies of two homozygous transgenic tomato lines, with different selection markers, expressing an intron-hairpin RNA cognate to the Tomato yellow leaf curl virus (TYLCV) C1 gene. The tested F6 and F4 progenies of the respective kanamycin- and basta-resistant plants exhibited unchanged field resistance to TYLCV and stably expressed the transgene-derived short interfering RNA (siRNAs) to represent 6 to 8% of the total plant small RNAs. This value outnumbered the average percentage of viral siRNAs in the nontransformed plants exposed to TYLCV-infested whiteflies. As a result of the RNAi transgene expression, a common set of up- and downregulated genes was revealed in the transcriptome profile of the plants selected from either of the two transgenic events. A previously unidentified geminivirus causing no symptoms of viral disease was detected in some of the transgenic plants. The novel virus acquired V1 and V2 genes from TYLCV and C1, C2, C3, and C4 genes from a distantly related geminivirus and, thereby, it could evade the repressive sequence-specific action of transgene-derived siRNAs. Our findings shed light on the mechanisms of siRNA-directed antiviral silencing in transgenic plants and highlight the applicability limitations of this technology as it may alter the transcriptional pattern of nontarget genes.

Embryonic gene expression of Coregonus palaea (whitefish) under pathogen stress as analyzed by high-throughput RNA-sequencing.

Most fishes produce free-living embryos that are exposed to environmental stressors immediately following fertilization, including pathogenic microorganisms. Initial immune protection of embryos involves the chorion, as a protective barrier, and maternally-allocated antimicrobial compounds. At later developmental stages, host-genetic effects influence susceptibility and tolerance, suggesting a direct interaction between embryo genes and pathogens. So far, only a few host genes could be identified that correlate with embryonic survival under pathogen stress in salmonids. Here, we utilized high-throughput RNA-sequencing in order to describe the transcriptional response of a non-model fish, the Alpine whitefish Coregonus palaea, to infection, both in terms of host genes that are likely manipulated by the pathogen, and those involved in an early putative immune response. Embryos were produced in vitro, raised individually, and exposed at the late-eyed stage to a virulent strain of the opportunistic fish pathogen Pseudomonas fluorescens. The pseudomonad increased embryonic mortality and affected gene expression substantially. For example, essential, upregulated metabolic pathways in embryos under pathogen stress included ion binding pathways, aminoacyl-tRNA-biosynthesis, and the production of arginine and proline, most probably mediated by the pathogen for its proliferation. Most prominently downregulated transcripts comprised the biosynthesis of unsaturated fatty acids, the citrate cycle, and various isoforms of b-cell transcription factors. These factors have been shown to play a significant role in host blood cell differentiation and renewal. With regard to specific immune functions, differentially expressed transcripts mapped to the complement cascade, MHC class I and II, TNF-alpha, and T-cell differentiation proteins. The results of this study reveal insights into how P. fluorescens impairs the development of whitefish embryos and set a foundation for future studies investigating host pathogen interactions in fish embryos.

Genome analyses suggest the presence of polyploidy and recent human-driven expansions in eight global populations of the honeybee pathogen Nosema ceranae.

Nosema ceranae is a microsporidian pathogen whose infections have been associated with recent global declines in the populations of western honeybees (Apis mellifera). Despite the outstanding economic and ecological threat that N. ceranae may represent for honeybees worldwide, many aspects of its biology, including its mode of reproduction, propagation and ploidy, are either very unclear or unknown. In the present study, we set to gain knowledge in these biological aspects by re-sequencing the genome of eight isolates (i.e. a population of spores isolated from one single beehive) of this species harvested from eight geographically distant beehives, and by investigating their level of polymorphism. Consistent with previous analyses performed using single gene sequences, our analyses uncovered the presence of very high genetic diversity within each isolate, but also very little hive-specific polymorphism. Surprisingly, the nature, location and distribution of this genetic variation suggest that beehives around the globe are infected by a population of N. ceranae cells that may be polyploid (4n or more), and possibly clonal. Lastly, phylogenetic analyses based on genome-wide single-nucleotide polymorphism data extracted from these parasites and mitochondrial sequences from their hosts all failed to support the current geographical structure of our isolates.

Sequencing and characterizing the genome of Estrella lausannensis as an undergraduate project: training students and biological insights.

With the widespread availability of high-throughput sequencing technologies, sequencing projects have become pervasive in the molecular life sciences. The huge bulk of data generated daily must be analyzed further by biologists with skills in bioinformatics and by "embedded bioinformaticians," i.e., bioinformaticians integrated in wet lab research groups. Thus, students interested in molecular life sciences must be trained in the main steps of genomics: sequencing, assembly, annotation and analysis. To reach that goal, a practical course has been set up for master students at the University of Lausanne: the "Sequence a genome" class. At the beginning of the academic year, a few bacterial species whose genome is unknown are provided to the students, who sequence and assemble the genome(s) and perform manual annotation. Here, we report the progress of the first class from September 2010 to June 2011 and the results obtained by seven master students who specifically assembled and annotated the genome of Estrella lausannensis, an obligate intracellular bacterium related to Chlamydia. The draft genome of Estrella is composed of 29 scaffolds encompassing 2,819,825 bp that encode for 2233 putative proteins. Estrella also possesses a 9136 bp plasmid that encodes for 14 genes, among which we found an integrase and a toxin/antitoxin module. Like all other members of the Chlamydiales order, Estrella possesses a highly conserved type III secretion system, considered as a key virulence factor. The annotation of the Estrella genome also allowed the characterization of the metabolic abilities of this strictly intracellular bacterium. Altogether, the students provided the scientific community with the Estrella genome sequence and a preliminary understanding of the biology of this recently-discovered bacterial genus, while learning to use cutting-edge technologies for sequencing and to perform bioinformatics analyses.

bacterial infection drives the expression dynamics of microRNAs and their isomiRs.

The optimal coordination of the transcriptional response of host cells to infection is essential for establishing appropriate immunological outcomes. In this context, the role of microRNAs (miRNAs)--important epigenetic regulators of gene expression--in regulating mammalian immune systems is increasingly well recognised. However, the expression dynamics of miRNAs, and that of their isoforms, in response to infection remains largely unexplored. Here, we characterized the genome-wide miRNA transcriptional responses of human dendritic cells, over time, to various mycobacteria differing in their virulence as well as to other bacteria outside the genus Mycobacterium, using small RNA-sequencing. We detected the presence of a core temporal response to infection, shared across bacteria, comprising 49 miRNAs, highlighting a set of miRNAs that may play an essential role in the regulation of basic cellular responses to stress. Despite such broadly shared expression dynamics, we identified specific elements of variation in the miRNA response to infection across bacteria, including a virulence-dependent induction of the miR-132/212 family in response to mycobacterial infections. We also found that infection has a strong impact on both the relative abundance of the miRNA hairpin arms and the expression dynamics of miRNA isoforms. That we observed broadly consistent changes in relative arm expression and isomiR distribution across bacteria suggests that this additional, internal layer of variability in miRNA responses represents an additional source of subtle miRNA-mediated regulation upon infection. Collectively, this study increases our understanding of the dynamism and role of miRNAs in response to bacterial infection, revealing novel features of their internal variability and identifying candidate miRNAs that may contribute to differences in the pathogenicity of mycobacterial infections.